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1.  DERIVATION  OF  TARGET  POSITION  EQUATIONS 


Before  an  analysis  of  target  position  error  can  be  performed,  the  equa- 
tions defining  the  problem  geometry  must  be  derived.  In  the  following  deriva- 
tion, the  target  position  is  described  as  a function  of  those  parameters  which 
are  both  pertinent  to  the  problem  geometry  and  capable  of  being  measured  by  on- 
board equipment  in  the  Scout  helicopter.  The  following  parameters  are  con- 
sidered: Scout  UTM  coordinates.  Scout  altitude.  Scout  heading.  Scout  attitude 

with  respect  to  local  vertical,  azimuth  and  elevation  angles  of  the  target 
relative  to  the  Scout  body  axis,  and  slant  range  from  the  Scout  to  the  target. 


S 


(S). 


Let  the  axes  Xq,  Yq,  and  ZQ 
From  Figure  1,  we  have 

Xqt  = Rfs  • cos  e . cos  a 

Yqt  = Rts  • cos  e • s*n  a 
zqt  - Rts  • sin  e 


represent  the  body  axes  of  the  Scout  aircraft 

la 

lb  (1) 

lc 


where  Xqx»  Yqt,  and  Zqt  are  the  target  (T)  coordinates  in  the  XYZ  system,  a 
and  e are  the  relative  azimuth  and  elevation,  respectively,  of  the  target  with 
respect  to  the  Scout's  body  axis  and  Rts  Is  the  slant  range  between  the  Scout 
and  target.  Consider  now  the  Scout's  attitude  with  respect  to  a horizontal 
reference  system  X3,  Y^,  Z3  with  X3  pointing  North,  Y3  pointing  East,  and  Z3 
pointing  vertically  down.  Take  first,  a rotation  through  an  angle  H (heading) 
about  the  Z3  axis  as  illustrated  in  Figure  2. 
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Figure  2.  Heading  angle  rotation 


Then 


— 1 
X 

l 

cos  H -sin  H 0 

"x2- 

*3 

= 

sin  H cos  H 0 

_Z3_ 

0 0 1 

z2 

(2) 


Second,  perform  a rotation  through  an  angle  E (elevation)  about  the  Y2  axis, 
as  illustrated  in  Figure  3. 


then 


x2 

Y2 

_ 

cos  E 0 

0 1 

sin  E 
0 

1 

H iH 

X 5* 

1 

Z2 

-sin  E 0 

cos  E 

Z1 

(3) 


Finally,  a rotation  about  an  angle  <}>  (roll)  is  performed  about  the  axis  as 
illustrated  in  Figure  4. 


X, , X, 


Figure  4.  Roll  angle  rotation 


Then 

~1  0 0 
0 cos  <f>  -sin  <J> 

0 sin  <{)  cos  <J> 

Combining  equations  2,  3,  and  4: 


~xl~ 

Y1 

= 

Zl 

x0 


X3T 

CH 

-SH 

0 

_CE 

0 

SE 

1 

0 

o’" 

Y3T 

as 

SH 

CH 

0 

0 

1 

0 

0 

C<t> 

-Sep 

Z3T 

0 

0 

1 

-SE 

0 

CE 

0 

S(J) 

C(j) 

where  CH 

denotes 

cos  (H) , 

etc. 

, , for  brevity. 

Multiplying 

the 

Xqt 

Y 

OT 

50T 


(4) 


(5) 


substituting  equations  la,  lb,  and  lc  into  (5)  yields 
CHCE 


1 

X 

- U 

y3T 

= 

Z3T 

SHCE 


-SE 


CHSES<J> 

-SHC$ 

SHSES<f> 

-CHC<f> 

CES4> 


CHSEC<}) 

-f-SHStf) 

SHSEC4) 

-CHSp 

CEC<() 


Rts  Ceca 

R'j'S  CEScx 

Ris  Se 

(6) 
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y (NORTH) 


Figure  5.  Definition  of  ground  coordinate  system 

If  the  X3,  Y3,  Z3  coordinate  system  is  now  referenced  to  a horizontal 
ground  coordinate  system,  X,  Y,  Z with  X East,  Y North,  and  Z up,  as  shown  in 
Figure  5 there  results 

XT  = Xs  + Y3T  7a 

YT  = Ys  + X3T  7b  (7) 

Z^  — Zg  — Z3T  7c 


where  Xx>  Yx,  Zx  are  the  target  coordinates  and  Xg,  Yg,  Zg  are  the  Scout  co- 
ordinates in  the  X,  Y,  Z system.  Combining  equations  (7a),  (7b),  and  (7c)  with 
equation  (6)  yields 


"yt" 

“Ys" 

CHCE 

+CHSES<f> 

-SHC<j> 

+CHSEC  <p 
+SHS<p 

CeCa 

XT. 

= 

Xs 

+ rts 

SHOE 

+SHSES<f> 

+CHC<}> 

+SHSEC<|> 

—CHS 

CeSa 

— Z»p 

-ZS 

-SE 

CES<f> 

CEC<(> 

Se 

Assuming  that  X and  Y represent  UTM  coordinates  and  Z represents  altitude, 
we  now  have  an  equation  defining  the  target  location  as  a function  of 

9 Scout  UTM  coordinates  Xs  and  Ys 

• Scout  altitude  Zs 


4 


Scout  heading,  H 


• Scout  attitude  with  respect  to  local  vertical  (elevation,  E and 

roll,  <J>) 

• Azimuth  angle,  a,  of  the  target  relative  to  the  Scout 

• Elevation  angle,  e,  of  the  target  relative  to  the  Scout 
m Slant  range  , RTS  from  the  Scout  to  the  target 

i.e..  Target  Location  * f (Xs,  Yx,  Zs,  H,  E,  <j),  a,  e,  Rxs) 

In  component  form,  equation  (8)  becomes 

XT  “ xs  + RTS  fSH  [Ce(CECa  + SES4>Sa)  + SEC<J>Se] 

+ CH(C4>CeSa  - s<j>Se)} 

yT  « Ys  + Rxs  CH{[Ce(CECa  + SES<j>Sa)  + SEC4>Se] 

- SH(C<f>CeSa  - S4.Se)} 

ZT  - zst+rts  {Ce(SECa  - CES4>Sa)  - CEC4>Se}  9c 


9a 


9b 


(9) 


It  is  worth  nothing  here  that  equation  (8)  may  be  rewritten  as  follows: 


1 

RTs 


H 

CO 

1 

CeCa 

xT  - xs 

= M 

CeSa 

-(zT  - Zs)_ 

se 

, M = matrix  in  equation  (8) 


or 


1 

rTS 


"yt  - 

Ys" 

CeCa 

mT 

XT  - 

xs 

s 

CeCa 

(Z-p  - 

zs)_ 

Se 

(8a) 


Then,  if  we  pick  X 
can  be  calculated  as  RTS 
e can  be  determined  by  equation  (8a5  above. 


, Y , Z , XT,  Yt,  Z , H,  E.  and  4>,  the  variable  Rts 
| = I(XT  - X )2  + (Y^  - Ys)2  + (zT  - Zs)2]l/2  and  a and 


2.  DERIVATION  OF  ERROR  EQUATIONS 


A linear  error  analysis  was  performed  on  the  target  position  variables 
XT,  Y^,  and  Zj.  Considering,  first,  Xx,  we  write 

Xt  = Xx  (Xg , Yg,  Zg , H,  E,  <j> , a,  £,  Rxs) 

letting 

h = Xs 
C2  * Ys 

^3  " Zs 
U = H 

^5  " E 
?6  = + 
h = a 

58  - e 
% = Rts 


and  taking  the  total  derivative  of  X.j.  yields 


3XT  dfcL 

3^ 


we  now  define  dx^  = Ex,  the  error  in  the  computed  value  of  the  target  coordin- 
ate and  d£i  = e^,  the  error  in  the  measurement  of  the  Ci^h  variable. 


Then 


Ex  = Z ^T  e± 


Similarly 

9 a 

Ey  = Z fYr  e± 

i=l  3? 


Ez 


9 U* 
i=l  ^ 


10a 


10b  (10) 


10c 


Calculate  first  the  partial  derivatives  of  Xx  with  respect  to  the  5-j_ : 


_ 
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1.  axT/axs  = 1 

2.  3xx/3Ys  = o 

3.  3XX/3ZS  = 0 

4.  3XX/3H  = Rts  {CH[Ce(CECa +SES(})Sa)  +SEC<j)Se]  -SH(C<J)CeSa  - S<J>Se)} 


but,  inspection  of  (9b)  shows  this  is  just  Yx  -Yg 


3XX 

3H~ 


Yx  -Ys 


5.  3XX/3E  = Rxs(SH[Ce  (SECa  -CES<J>Sa)  -CEC<J>Se]} 
but,  inspection  of  (9c)  shows  this  is  just  SH(ZX  -Zg) 

3XX/3E  = SH(ZX  -Zs) 

6.  3XT/B<p  = RXS(SH  (CeSEC<J>Sa  -SES4>S  e]  +CH[-S«{)CeSa  — C4>S  e]  } 

7.  3XX/3a  = RXS(SH[CE(-CESa  +SES<t>Ca)  ] +CHC4>CeCct} 

8.  3Xx/3e  = Rxs{SH[-Se(CECa  +SES<t>Sa)  +SEC<J>Ce]  +CH(-C$SeSa  -S<J>Ce)  } 

9.  3XX/3RTs  = (SH[  (CECa  +SES<{>Sa)Ce  +SEC<{>Se]  +CH[C<|>CeSa  -S<J>Se]} 

but,  inspection  of  (9a)  shows  this  is  just  (XX  -Xg)/Rxs 

1*1  = (Xx  -XS)/RTS 
3 Rts 

Next,  calculate  the  partial  derivatives  of  Yx  with  respect  to  the  £i: 

1.  3YX/3X$  = 0 

2.  3YX/3YS  = 1 

3.  3YX/3ZS  = 0 

4.  3YX/3E  = - (Xx  -Xs) 

5.  3YX/9E  = CH(ZX  -Zs) 

6.  3YX/ 3 = RXg{CH  [CeSEC^Sa  -SES$>Se]  -SH[-S<frCeSa  -C<{>Se]} 

7.  3YX/3a  = Rxs{CH[Ce(-CESa  +SES<t>Ca)]  -SH(C4>C eCa)  } 

8.  3Yx/3e  = RXg{CH[-S e(CECo  +SES<f>Sct)  +SEC<f>Ce]  -SH[-C<{>SeSa  -S<J>Ce]} 

9.  3YX/3RXS  - (YX  -Ys)/RXS 
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Finally,  calculate  the  partial  derivatives  of  Zj  with  respect  to  the 


1.  3ZT/3XS  = 0 

2.  3Zx/3Yg  = 0 

3.  3ZX/3ZS  = 1 

4.  3ZT/3H  = 0 

5.  3Zt/3E  = -SH(Xt  -Xs)+CH(Yt  -Ys) 

6.  3zt/34>  = RTS(Ce(CEC(j)Sa)  -CES<j>Se} 

7.  3ZT/3a  = -Rts  Ce(CES<j>Ca  +SESa) 

8.  3Zx/3e  = -RTs(-Se(CES4)Sa  -SECa)  +CEC<f>Ce} 

9.  3ZT^RtS=  (ZT  - Zs)/RTS 


Letting 

Xts  = Yx  -Xs 
YTS  = Yt  -Yg 

ZTS  = ZT  ~zs 


There  results 


Ex 

1 

0 

0 

yts 

(CH)Zts 

ey 

= 

0 

1 

0 

-xTS 

(SH)Zts 

EZ 

0 

0 

1 

0 

— (SH)Xt’c 
-(CH)YXS 

3X,. 

9xt 

3Xt 

XTS 

W 

3a 

3eT 

rts 

3YX 

8Y*p 

3Y»p 

yts 

3<J) 

3a 

rXs 

3ZX 

3ZX 

3Z«p 

ZTS 

w 

3cT 

3e 

eE. 

e<P 


eR 

(ID 
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The  matrix  equation  11  defines  the  total  error  in  target  position  as  a 
function  of  the  nine  Scout-target  variables  Xs,  Ys,  Zs,  H,  E,  <j>,  a,  e,  and 
R^S  and  the  nine  associated  errors,  e^. 

3.  ANALYSIS  OF  THE  ERROR  EQUATIONS 

Assume  that  the  nine  Scout- target  variables  have  been  assigned  specific 
values  and,  therefore,  the  matrix  elements  in  equation  (11)  have  been  fixed. 
Let  these  elements  be  denoted  C..  (i  = row;  j = column) . Next,  assume  that 
the  ei  are  independent,  normal lyJ distributed  random  variables  with  density 
functions. 

-(n  -y1)2/2ai2 

fi(n)  = N[Ui,  Oil  = e 

a ±v2tT 


It  can  then  be  shown  that  the  density  functions  corresponding  to  the  random 
variables 

9 

Ex  = ^ Cy  ej 

9 

Ey  * E ^21 

j=l 

9 

Ez  = ^ Co.«  ®1 

j-1  2 

are  normal  with  means 


k = 1,  2,  3 


and  standard  deviations 


2\l/2 


k - 1,  2,  3 


That  is,  the  normal  density  functions  corresponding  to  Ey,  EY,  and  Ez  are 


9 


um 


— 


9 

9 

fE*  = N 

AClj  ^ 

•.z.(cij  aiy 
J-l 

9 

9 

% " N 

E Co-5  ^4  * 
j-1  J J 

9 

9 

55 

II 

N 

W 

4-1 

E Cjj  \1*  , 

j-i  L 

£ (Coj  0.1) 
j-1  J J 

1/2 

• 1/2 

J 

1/2 


Suppose,  now,  that  it  is  desired  to  perform  the  above  analysis  for  many 
Scout-target  states.  For  simplicity,  only  the  distributions  in  Ex  will  be  con- 
sidered. For  N Scout-target  states,  there  results  N normally  independent 
probability  distribution  functions  (pdf)  f^,  with  corresponding  means  y^;  and 
standard  deviations  a^,  i - 1,  2,  ...,  N.  We  now  pose  the  following  question: 


If  a sample  is  taken  at  random  (with  equal  likelihood)  from  one  of  the 
distributions,  what  is  the  distribution  of  the  sample? 


The  sample  space  is  shown  in  Figure  6. 


Figure  6.  N normally  independent  probability  distribution 
functions  comprising  sample  space 


Since  the  sample  is  taken  with  equal  likelihood,  the  sample  space  is  com- 
posed of  N functions  of  the  form: 


fz(i,y) 


-(y  ~yi)2 

2°i 


V^TT  Oj[  N i = 1,  2 N 

Summing  over  i yields  the  pdf  of  the  sample: 

-(y  -Pi)2 


N 

My)  - i i 

/2tt  N i=l 


2 0±2 


Then  the  mean  of  the  sample  may  be  calculated  as 

00 

Uz  = f Zfz(Z)dZ 

-TOO 

(z  -yjL)2/2ai2 


N 


~ 1 Z 1 f Ze 

N i=l  m 0±  -L 

N 

= i 2 y, 

N ±ml  1 


dz 


and  the  variance  as 


00 

oz2  =J  z2  fz(z)dz  - yz2 


N 


i Z _1_  f 
N i=1v^f  0±  J°° 


Z2e 


~(z  -Ui)2/2a:L2 


dz  "^z 


<V 


I £ fa,*  <V>  -11, 


4.  RESULTS 

Several  techniques  were  used  in  performing  the  analysis.  The  first  con- 
sisted of  choosing  nominal  values  for  Xs,  Ys,  Zj,  and  H and  allowing  X-p.  Y-j>> 
Zs,  E,  and'  <J>.  to  vary  in  increments  over  predetermined  ranges.  Rxs>  and  e 
were  calculated  according  to  equation  (8a).  For  convenience,  Xj  and  Yj  were 
transformed  into  polar  coordinates  by  the  equations: 

XT  = r.sin(e) 

Yt  = r.cos(6) 
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and  then  the  parameters  r and  0 were  varied.  The  following  tables  list  the 
standard  deviations  corresponding  to  each  e^  and  the  maximum,  minimum,  and  in- 
cremental values  of  the  Scout-target  parameters. 


TABLE  1.  STANDARD  DEVIATIONS  AND  MEANS  CORRESPONDING  TO  EACH  VARIABLE,  e± 


error 

a 

y 

ex 

55.0  m 

0 

ey 

55.0  m 

0 

ez 

2.0° 

0 

eH 

1.0° 

0 

*E 

0.5° 

0 

e4> 

0.5° 

0 

ea 

0.2  mrad 

0 

ee 

0.2  mrad 

0 

eR 

10.  0 m 

0 

TABLE  2.  MAXIMUM,  MINIMUM,  AND  INCREMENTAL  VALUES 
OF  THE  SCOUT-TARGET  PARAMETERS 


PARAMETER 

MIN 

MAX 

INC 

xs 

0 

0 

0 

Ys 

0 

0 

0 

zs 

10  m 

35.  m 

5.  m 

0 = tan-1(XT/YT) 

0° 

O 

O 

18° 

r = (XT2  + YT2)l/2 

2000  m 

4000  m 

400  m 

Zx 

0 

0 

0 

H 

45° 

45° 

0 

E 

5° 

10° 

1° 

<P 

-4° 

1° 

1° 

a 

-45.89° 

45.66° 

— 

e 

0.84° 

11.72° 

— 

Rts 

2000.02  m 

4000.15  m 

— 

Allowing  five  parameters  to  assume  6 values  each  yields  6^  combinations 
or  7776  Scout-target  states. 
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The  resulting  standard  deviation  about  the  target  location  were 


0F  = 67,26  meters 
aEy  = 67.26  meters 
0£z  * 26.87  meters 

and  the  Circular  Probable  Error  (CEP)  was 


CEP  = 1. 1774aEx  = 79.19  meters 

In  order  to  determine  the  contribution  of  individual  errors  to  the  total 
error,  nine  additional  runs  were  made  on  the  computer,  each  with  all  but  one 
of  the  nine  standard  deviations  zeroed.  That  is,  the  first  run  was  made  with 
ex  nonzero,  the  second  with  ey  nonzero,  etc. 

Letting  aEx^(e^)  * variance  of  target  x position  with  the  i t^1  error  non- 
zero, etc. 

2 

and  0EX  = variance  of  target  X position  with  all  errors  nonzero, 

etc. 


Then 


Vi (Ex) 


OEx^ 

°EX2 


measure  of  the  contribution  of  the  it^  error  to  the 
total  error  in  X. 


The  following  table  lists  the  results  of  the  above  calculations: 
TABLE  3.  CONTRIBUTION  OF  INDIVIDUAL  ERRORS  TO  THE  TOTAL  ERROR 


Nonzero 

Input  Error  (i) 

Vi (Ex) 
(%) 

Vi(Ey) 

(%) 

V±(EZ) 

(%) 

ex 

66.87 

0 

0 

ey 

0 

66.87 

0 

ez 

0 

0 

0.55 

eH 

31.87 

31.87 

0 

eE 

4.97  by  1CT4 

4.97  by  10-4 

75.52 

e4> 

0.16 

0.16 

23.88 

ea 

4.08  by  10-3 

4.08  by  10-3 

3.46  by  10-4 

e£ 

3.54  by  10"  5 

1.99  by  10-5 

0.05 

eRTS 

1.10 

1.10 

1.12  by  10~3 
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It  is  clear  from  Table  3 that,  within  the  range  of  Scout- target  state 
parameters  chosen,  the  navigation  errors  (ex  and  ey)  and  the  heading  error 
(ey)  are  the  major  contributors  to  the  target  location  error  in  the  x-y  plane, 
which  is  the  plane  of  the  UTM  grid.  For  this  reason  it  was  decided  to  take  a 
more  detailed  look  at  the  relationship  between  navigation  error,  heading  er- 
ror, and  total  error.  Realizing  that  the  contribution  to  the  total  error  from 
heading  error  is  dependent  on  range,  six  runs  were  made  with  all  parameters 
and  errors  as  in  Tables  1 and  2 except  that  the  parameter  r was  held  constant 
throughout  each  run.  The  results  of  each  run  were  then  a function  of  r which 
assumed  values  ranging  from  1,000  meters  for  the  first  run  to  6,000  meters  for 
the  sixth  run.  Table  4 shows  the  results  of  those  runs. 

TABLE  4.  CONTRIBUTION  OF  NAVIGATION  AND  HEADING  ERRORS  TO  TOTAL  ERROR 


r 

meters 

TOTAL 
ERROR  (m) 
(CTEX  = CTEY) 

ERROI 

NAVI 

l DUE  TO 
[GATION 

ERROR  DUE  TO 
HEADING 

ABSOLUTE 

(meters) 

RATIO  OF 
VARIANCES 

ABSOLUTE 

(meters) 

RATIO  OF 
VARIANCES 

1,000 

56.82 

55.0 

93.69% 

12.34 

4.72% 

2,000 

60.72 

55.0 

82.04% 

24.68 

16.52% 

3,000 

66.73 

55.0 

67.93% 

37.02 

30.78% 

4,000 

74.32 

55.0 

54.77% 

49.37 

44.13% 

5,000 

83.07 

55.0 

43.84% 

61.71 

55.18% 

6,000 

92.65 

55.0 

35.24% 

74.05 

63.88% 

Finally,  120  additional  runs  were  made  on  the  computer  each  with  a fixed 
r,  navigation  error,  and  heading  error.  The  intent  of  these  runs  was  to  show 
the  trade-off  between  heading  error  and  navigation  error  for  varying  ranges. 
The  results  are  listed  in  Tables  5 through  10. 


TABLE  5.  TRADE-OFF  BETWEEN  HEADING  ERROR  AND  NAVIGATION 
ERROR  FOR  RANGE  EQUAL  TO  1,000  METERS 


Navigation 

Error 

CEP  (meters) 

Resulting  Tgt  Loc  CEP 

Heading  Errors  (°) 

0.5 

1.0 

1.5 

2.0 

1,000 

1000.02 

1000.11 

1000.24 

1000.43 

500 

500.16 

500.32 

500.58 

500.96 

250 

250.21 

250.53 

251.04 

251.79 

100 

100.57 

101.36 

102.66 

104.43 

50 

51.25 

52.78 

55.22 

58.47 
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TABLE  6.  TRADE-OFF  BETWEEN  HEADING  ERROR  AND  NAVIGATION 
ERROR  OR  RANGE  EQUAL  TO  2,000  METERS 


Navigation 

Error 

CEP  (meters) 

Resulting  Tgt  Loc  CEP 

Heading  Errors  (°) 

0.5 

1.0 

1.5 

2.0 

1,000 

1000.11 

1000.43 

1000.95 

1001.68 

500 

500.32 

500.96 

502.00 

503.48 

250 

250.52 

251.79 

253.88 

256.78 

100 

101.37 

104.46 

109.39 

115.95 

50 

52.81 

58.51 

66.92 

77.18 

TABLE  7.  TRADE-OFF  BETWEEN  HEADING  ERROR  AND  NAVIGATION 
ERROR  FOR  RANGE  EQUAL  TO  3,000  METERS 


Navigation 

Error 

CEP  (meters) 

Resulting  Tgt  Loc  CEP 

Heading  Errors  (°) 

0.5 

1.0 

1.5 

2.0 

1,000 

1000.25 

1000.95 

1002.14 

1003.80 

500 

500.60 

502.02 

504.37 

507.66 

250 

251.07 

253.89 

258.52 

264.88 

100 

102.69 

109.42 

119.78 

132.94 

50 

55.30 

66.96 

82.82 

100.91 

TABLE  8.  TRADE-OFF  BETWEEN  HEADING  ERROR  AND  NAVIGATION 
ERROR  FOR  RANGE  EQUAL  TO  4,000  METERS 


Navigation 

Error 

CEP  (meters) 

Resulting  Tgt  Loc  CEP 

Heading  Errors  (°) 

0.5 

1.0 

1.5 

2.0 

1,000 

1000.42 

1001.70 

1003.80 

1006.74 

500 

500.97 

503.49 

507.67 

513.46 

250 

251.82 

256.80 

264.89 

275.83 

100 

104.52 

115.0 

132.96 

153.59 

50 

58.60 

77.25 

100.95 

126.90 
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TABLE  9.  TRADE-OFF  BETWEEN  HEADING  ERROR  AND  NAVIGATION 
ERROR  FOR  RANGE  EQUAL  TO  5, 000  METERS 


Navigation 

Error 

CEP  (meters) 

Resulting  Tgt  Loc  CEP 

Heading  Errors  (°) 

0.5 

1.0 

1.5 

2.0 

1,000 

1000.67 

1002.65 

1005.93 

1010.52 

500 

501.45 

505.39 

511.87 

520.82 

250 

252.78 

260.49 

272.86 

289.29 

100 

106.80 

123.96 

148.20 

176.63 

50 

62.60 

88.75 

120.32 

153.99 

TABLE  10.  TRADE-OFF  BETWEEN  HEADING  ERROR  AND  NAVIGATION 
ERROR  FOR  RANGE  EQUAL  TO  6,000  METERS 


Navigation 

Error 

CEP  (meters) 

Resulting  Tgt  Loc  CEP 

Heading  Errors  (°) 

0.5 

1.0 

1.5 

2.0 

1,000 

1000.97 

1003.82 

1008.54 

1015.11 

500 

502.04 

507.69 

516.96 

529.68 

250 

253.94 

264.93 

282.29 

304.95 

100 

109.53 

133.03 

164.93 

201.25 

50 

67.16 

101.04 

140.40 

181.70 

5.  CONCLUSIONS 


Within  the  context  of  this  report,  it  is  clear  that  the  navigation  and 
heading  system  errors  have  the  most  significant  impact  upon  target  location 
error.  This  is  demonstrated  by  the  data  presented  in  Tables  3 and  4.  Tables 
5 through  10  demonstrate  that,  for  large  navigation  errors,  the  target  loca- 
tion CEP  is  relatively  insensitive  to  changes  in  heading  errors  and  changes 
in  range.  For  smaller  navigation  errors,  the  contribution  from  heading  errors 
becomes  increasingly  significant  for  increasing  range.  As  seen  in  Table  10, 
for  instance,  with  r = 6,000  meters  and  navigation  error  = 50  meters,  the  tar- 
get location  error  increases  from  67.16  meters  to  181.7  meters  (a  factor  of 
2.7)  as  the  heading  error  increases  from  0.5°  to  2°. 
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Finally,  it  should  be  emphasized  that  when  using  the  results  presented 
in  this  report,  one  must  be  mindful  of  the  fact  that  all  system  errors,  in- 
cluding navigation  and  heading,  have  been  assumed  to  be  independent  when,  in 
the  real  world,  this  may  not  be  the  case. 

6.  ACKNOWLEDGEMENT 


The  author  wishes  to  express  his  graditude  to  Dr.  R.  J.  Niemela  for  his 
valuable  assistance  during  the  performance  of  this  study. 


APPENDIX 


PROGRAM  LISTINGS  AND  SAMPLE  OUTPUT 


COMMON  XS,  YS.ZS. XT. YT.ZT.H, E, PHI , ALPHA, EPS, RTS,  PD<3,9)  , IXYZ 
DIMENSION  AME AN ( 9) 

DIMENSION  V IN  IT (9), VF  I N AL < 9 ) , V I NC < 9 ) , ERROR < 9 ) ,NI (9) 

DIMENSION  ERINIK9)  ,ERFIN<9>  »ERINC<9)  ,Ml  (9) 

EQUIVALENCE  (NI (1) ,N1) , (NI (2) ,N2> , ( N I (3),N3>, (NI (4) ,N4) , (NI (5) ,N5) 
l»(Nl(6),N6),<NI(7),N7)t(NI(8),N8),(NI(9),N9),(Ml(l),Ml),(MI(2)»M2) 
2»<Ml(3).H3)»(Ml(4)fM4),(Ml(5),M5),(Ml(6),M6),(Ml(7)»M7),(M-!(8),M8> 
3,  (Ml (9) ,M9) 

DOUrtLE  PRECISION  SMS  I GX , SMS  I GY , SMS  I GZ . SUMEMX , SUMEMY , SUMEMZ 
LOGICAL  CREG 

RE AD ( 3 , 5 ) IXYZ.NL  I NES, ITEST 
5 FORMATf 11, 213) 

RE  AD  ( 3 , 10  ) (AMEANU  ),  1=1,9) 

READ ( 3,  10 ) (ERINITC  I ) , 1=1,9) 

RE AD ( 3 , 10 ) (ERFIN( I ), 1=1,9) 

RE AD ( 3 , IB ) (ERINCt I ), 1=1,9) 

10  FORMAT(3F6.1,6<1PE10.3>) 

REAQ ( 3, 10 ) (VINIT(I),I=1,9) 

READ( 3, 10 ) (VFINAL( I), 1=1.9) 

REAO ( 3 , 10 ) ( V I NC ( I ), 1 = 1,9) 

WRITE(l.ll) 

11  FORMAT( 1X.45HSET  CREG(l)  FOR  POLAR  X-Y  TARGET  COORDINATES.) 

PAUSE 

NTEST=0 
PD( 1 » 1 ) =1 . 0 
PD( 1,2)*0.0 
PO( 1 , 3 ) =0 . 0 
PD ( 2 , 1 ) =0 . 0 
PO( 2, 2) *1 . 0 
PD ( 2, 3 ) =0 . 0 
PO(3,1)=0.0 
PO(3,2)=0.0 
PD( 3, 3 ) *1 . 0 
PO ( 3 , 4 ) =0 , 0 
DO  20  1=1.9 

I F ( ASS ( ER I NC ( I ) ) .GT.  1.0E-20)  GO  TO  12 
MI  ( I ) = 1 
GO  TO  13 

12  MI  < I ) = ABS<  <ERFIN< I )-ERINlT( I ) )/ERINC< I > ) + l,5 

13  IF<A8S<VINC<  I ) ) .GT.  1.0E-20)  GO  TO  15 
NI  ( I )=1 

GO  TO  20 

15  NI  ( I ) =ABS ( ( VF  I N AL ( I ) - V I N I T ( I ) ) / V I NC ( I ) ) *1 . 5 
20  CONTINUE 

C0N*3. 14159/180.0 
KOUNT  = 0 

DO  30  I El«l » Ml 

ERROR(1)=ERINIT(1)+FLOAT( IE1-1)*ERINC(1) 

DO  30  IE2=1,M2 

ERR0R<2>=ERIMT<2)+FL0AT< I E2- 1 ) *ER I NC < 2 > 

DO  30  ' I E3*l , M3 

ERR0R(3)=ERINIT(3)*FL0AT( I E 3- 1 ) *ER I NC  < 3 ) 
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DO  30  IE4*1,M4 

ERROR (4)s<ERINIT(4)+FL0 AT ( IE4-1)*ERINC<4) ) *CON 
00  30  IE5=1,M5 

ERROR(5)=(ERINIT(5)+FLOAT( IE5-1 )*ERINC(5) )*CON 
DO  30  I Eft  = 1 * M6 

ERR0R{6)s(ERINIT(6)*FL0AT( I E6-1 ) *6R INC ( 6 ) )*CON 
DO  30  IE7=1,M7 

ERROR(7)sERINIT(7)+FLOAT( I E 7 — 1 ) *ER I NC ( 7 ) 

DO  33  IE8*1,M8 

ERR0R(8)=ERINIT(8)+FL0AT( I E 8 — 1 ) *ER I NC ( 8 ) 

00  30  I E9  = 1 * M9 

ERROR (9 )=ER IN  IT (9) ♦FLOAT* I E9-1 ) *ER I NC ( 9 ) 

22  — NSTATE=0 
SUMEMX=0. 0 
SUMEMYs0. 0 
SUHEMZ=0. 0 
SMSIGXs0.0 
SMS  I GY  = 0 . 0 
SMSIGZ=0.0 
ALPMIN=1. 0E*6 
ALPMAX=0.0 
EPSMlN=1.0E+6 
EPSMAXs0, 0 
RNGM I N*1 . 0E*6 
RN&MAXr0.0 
DO  35  1*1, N1 

XS*V!NIT< 1) ♦FLOAT* I-1)*VINC<1> 

DO  35  J*1,N2 

YS»V INI T<2> ♦FLOAT* J-1)*VINC<2) 

DO  35  K*1,N3 

2S*VINIT<  3-)^FL0AT(K-l)*VlNC(3) 

DO  35  L*1,N4 

I NI T<4)^FtO AT (L-1)*VINC( 4 >)*CON 
HP=H/CON 
DO  35  M»1,N5 

Er( VINIT(5)+FL0AT(M-1)*VINC(5) )*CON 

tP^e/CON 

DO  35  N*1.N6 

PRt»iVIN-IT<6-)*FL0AT(N-l>*VlNC<6)  )*CON 

PHIP=PHI/C0N 

DO  35  I I*1,N7 

XTsVlNIT(7)*FL0AT( I I-1)*VINC<7> 
tFtlXYZ  .EQ.  0)  ALPHAsXT  *CON 
DO  35  JJ*1,N8 

-W  I-NITt8>^fL0AT(  J J-l)  *V  INC  ( 8) 

IF(CREC(1) ) X T= ( V I N I T ( 7 ) + FLOAT ( I I -1 ) *V I NC ( 7 ) ) *S I N { YT*CON ) 
I F ( CREG ( 1 ) ) YT* ( VINIT( 7) ♦FLOATt I 1-1 >*VINC(7) ) *COS ( YT*CON > 
I F ( I XYZ  .EQ.  0)  EPS*YT*CON 
DD  35  KK*1,N9 

ZT=VINIT(9)*FL0AT(KK-1)*VINC(9) 

I F ( I XYZ  .EQ.  0)  RTS=ZT 
NSTATEsNSTATE+1 
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CALL  GETPO 
ALPHAPsALPMA/CON 
EPSP=EPS/CON 
EX  = 0 . 0 

Ey«0.0 
EZ  = 0 . 0 
EMX=0. 

EMY=0. 

Emz= 0 . 

00  25  LL=1,9 

EX  = EX<*-P0(1,LL)*PD(1,LL)*ERR0R(LL)»ERR0R(LL) 

EY  s EY^PD(2,LL)*P0(2.LL)*ERR0R(LL)*ERR0R(LL) 

EZ  = EZ+PD(3,LL)*P0( 3,LL)*ERR0R(LL)*ERR0R(LL) 

EMX  = EMX*PD(1.LL)*AMEAN(LL) 

EMY  = EMY*PD<2,LL)*AMEAN<|_L> 

25  EMZ  = EMZ*P0(3,LL)*AMEAN(LL) 

SMSIGX=SMSIGX+EX 
SMS IGYsSMSI GY+EY 

smsigz=smsigz+ez 

SUMEMXaSUMEMX'fEMX 

SUMEMY=SUMEMY*EMY 

SUMEMZ»SUMEMZ+EMZ 

I F ( ALPM AP  .LT.  ALPMIN)  ALPM INsALPHAP 
IF( alphap.gt.  ALPMAX)  alpmax» alphap 
IF(EPSP  .LT.  EPSMIN)  EPSMIN=EPSP 
IF ( EPSP  .GT.  EPSMAX)  EPSM  AX  *EPSP 
IF (RTS  .LT.  RNGMIN)  RNGM I NsRTS 
IF (RTS  .GT.  RNGM AX ) RNGMAX=RTS 
35  CONTINUE 

STATE =SQRT!FLOAT< NS T4TE)) 
emx=sumemx/state 

EMY=SUMEMY/STATfe 
EMZ=SUMEMZ/STATb 
SIGXsOSQRT(5MSIGX)/STATE 
SIGY=0SQRT(SMSIGY)/STATE 
SIGZ  = OSQRT(SMSIGZ)/STATt: 

00  55  111=4,6 

55  ERROR!  U1)=ERR0R<  UD/CON 
WR I TE ( 6 » 100 ) 

100  FORMAT! 1H1.41X.44HI  NPUT  MATRIX  PARAMETERS/// 
C18X, 14HSC0UT  POSITION, 16X, 15HTARG6T  POS I T I ON , 5X , 7HHE AO  I NG , 5X , 5HP I T 
CCH.6X,  4HR0LL, 2X.8HREL . AZ  . . 2X.8HREL . EL  ., 5X , 5HRANGE/2 IX , 8H< METERS  ) 
C,  21 X,  8H  ( METERS ),llX,5H(0fcG) , 6X , 5H ( DEG ) , 5X , 5H ! OEG ) ,3X,5H(DEG) ,5X,5H 
C(DEG) ,9X, 3H!M)//13X,2HXS,8X.2HYS,8X,2HZS,8X,1HR,5X,5HTHETA,8X,2HZT 
C//) 

WRITE  (6, 110)  ( VINIT! I JK)  , I JK*1,3) , !VINIT( I JK) , I JK*7,9) , (VINIT! I JK) 
C,  I JK  = 4,6) , ALPMIN, EPSM IN, RNGM IN 
110  FORMAT! IX, 3HMIN, IX, 12F10.2) 

WRITE (6, 112)  (VFINAL! I JK) , I JK«l,3) , (VFINAL! I JK> , I JK*7,9) , tVFlNALC I 
CJK)  , I JK  = 4,6), ALPM AX, EPSMAX. RNGMAX 
112  FORMAT(1X,3HMAX,1X,12F10.2> 

WRITE (6, 114)  ( VINC! I JK) , I JK*1»3) , < VINC! I JK) , I JK*7,9) , (VING  ( I JK) , I 
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C JK  = 4 , 6 ) 

114  FORMAT(1X,3MINC,1X,9F10.2,3(6X,4H****)  ) 

WR I TE ! 6 » 116 ) NSTATE 

116  F0RMAK/1X,  19HNUMPER  OF  STATES  * ,17///) 

WR I TE ( 6 » 120 ) 

120  FORMAT ! 46X , 41H I NPUT  ERROR  STATISTIC  S////32X.16H 
CSCOUT  NAVIGATION, 4X»3HALT,3X,7HHEADING,5X,5HPITCH,6X,4HR0LL*2X,8HR 
CEL.  A2.  , 2X.8HREL.  EL ., 5X , 5HRANGE/36X , 8H ! METERS ), 8x , 3H ! M ) ,4X,5H 

C(DEG) ,6X,5H( DEG) ,5X,5H! DEG) ,3X,5H( RAD) ,5X,5H( RAD) ,4X,8H( METERS ) ,/3 
C4X,1MX,9X,1HY//) 

WR I TE ( 6 « 130  ) ( AMEAN( I JK) , I JK  = 1, 9)  , ( E-RROR  ( I JK  ) , I JK=  1 , 9 ) 

130  FORMAT (6X.4HME AN, 15X, 6F1 0. 1, 2F 10.4, F10.1/6X.18H3TANDARO  OEVlATION, 
C1X,6F10.1,2F10.4,F10.1/////) 

WRITE (6, 140) 

140  FORMAT! 16 X, 88 HCOMPU TED  ERROR  STATISTICS  FOR  RESULTING  DENSITY  FUNC 
1TI0NS  ABOUT  TARGET  X,Y,Z  COORD  I NATES//38X , 7HX  ERR0R.3X.7HY  ERROR, 3 
2X.7HZ  ERROR//) 

WRI TE (6, 150)  EMX.EMY, EMZ. SIGX.SIGY, SIGZ 

150  FORMAT! 15 X, 4HME AN, 16X, 3F10. 2/15 X, 18HST  AND ARD  OE V I AT  I ON , 2X , 3Fl 0. 2/ 1 
1 HI ) 

DO  155  122=4,6 

155  ERROR!  I22)=ERR0R( I22>*C0N 

30  CONTINUE 
CALL  EXIT 
ENO 
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SUBROUTINE  GETPD 

COMMON  XS»  YS.ZS#XT,YTf ZT.H,E.PHI,ALPHA»EPS#RTS»PD(3»9>,  IXYZ 
SH=SIN(H) 

CH=COS(H) 

SEx-SIN(E) 

CE=COS ( E ) 

SPHI=SINCPHI ) 

CPH  I = COS (PHI) 

IF ( IXYZ  .EQ.  0)  GO  TO  40 

XTS=XT-XS 

YTS=YT-YS 

ZTS=ZT-ZS 

RTS=SQRT(XTS*XTS+YTS*YTS+ZTS»ZTS) 

Tl=(-CH*SE*CPHI+SH*SPHI )*YTS 
T2=(SH*SE*CPHI+CH#SPHI )«XTS 
SEPS=(T1-T2-CE*CPHI*ZTS)/RTS 

CEPS2SQRT(1,0-SEPS«SEPS)  - - 

I F ( ABS ( CEPS ) .GT.  1.0E-10)  GO  TO  10 
EPS=1. 570796 
GO  TO  20 

10  EPS  = ATAN2(SEPS,CEPS)  - “ - — 

20  Tls-( CH*SE*SPHI ♦SH*CPHI )*YTS 

T2»CSH*SE*SPHI-CH«»CPW1-**XTS~  ' 

SALPHAs(Tl-T2-CE*SPHI*ZTS)/(RTS*CEPS) 

CALPHAsSORT(t.0-SALPHA^^ALPHAy  — 

IF( ABS(CALPHA)  .GT.  1.0E-10)  GO  TO  30 

ALPHA=1. 570796  — 

GO  TO  50 

30  ALPHA  = ATA-N2<SAIPHA,CALPHA-) - : — 

GO  TO  50 

4tf  SEPS=SINf€PS>  — ' : 

CEPS=COS ( EPS ) 

S AL.P*A*S  TN  t ALPHA  

CALPHA=COS( ALPHA) 

50  TlsCE*C*LPHA--Se*SPHl’»SAtPHA 
T2=SE*CPHI*SEPS 

T3xCE«SPHI  «SAtPHA+SE»CAtPK*r 

T4=CE#CPHI*SEPS 

T5T:CPHI*CEPS*SAL-P+T4-SPHlirSEPS 

T6=-SPHI*CEPS*SALPHA-CPHI*SEPS 

T7=-CE«SALPHA-SE*SPHI*CALPHA  - ■ 

T8=CPHI*CEPS*CALPHA 

T 9=CE«  C~A1_P  KA  p*  SPH  I*^A1_PWA 

T10=SE*CPHI*CEPS 

Tlls^CPHi*SEPSvSALPHA-SPHl*CtPS 

T12=CE*SPHI*CALPHA-SE*SALPHA 

T13*CE*CPH1*CEPS  - - ~ 

IF(IXYZ  .EQ.  1)  GO  TO  60 

XTxXS*RTS*(SH*(CEPS*T1-T2>*CH*T5>- 

YT=YS+RTS*(CH*(CEPS*Tl-T2)-SH*T5) 

ZTsZS-RTS*<C5PS*T3*T4>  

XTS=XT-XS 
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YTS-YT-YS 
ZTSsZT-ZS 
60  PD(1,4)=YTS 

PD ( 1#  5 ) *SH*ZTS 

PD<1»6)®RTS*(SH*(-SE*T5)+CH*T6) 

PD( 1, 7)=RTS*(SH*CPPS*T7+CM*T8) 

PD(1#  8)sRTS*(SH*( -SEPS*T9-T10>*CH*T11 > 
PD ( 1»  9 ) *XTS/RTS 
PD ( 2 » 4 ) *-XTS 
PD ( 2*  5 ) =CH*ZTS 

PD(2, 6)=RTS*(CH*(-SE*T5)-SH*T6) 
PD(2,7)sRTS*(CH*(CEPS»T7)-sH*T8) 
PD<2,8)sRTS*(CH*(-SEPS*T9-T10)-SH*T11 ) 
PD<  2,9)=YTS/RTS 
P0(3,5>s-(SH*XTS*CH*YTS) 
PD(3,6)=-RTS*CE*T5 
PD(3,7)*-RTS*CEPS*T12 
PD (3, 8)=-RTS*(-SEPS*T3+T13) 

PD( 3»9>»ZTS/RTS 

RETURN 

END 
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91-»T80I-WjI-VSIH 


INPUT  MATR  IX 


P A R A M E T £ R S 


N3 

4S 


SCOUT  POSITION 
(METERS  > 

TARGET  POSITION 
( MET EPS ) 

HEADING 

(OEG) 

PITCH 

(DEG) 

roll 

(OEG) 

REL.  az. 
(DEG) 

REL,  EL. 
(OEG) 

RANGE 

(M) 

~Y5 

TS 

75 

R 

THTT  a 2T 

MJN 

MAY 

INC 

0,00 

0,00 

0,00 

0.00 

0,00 

0,00 

10,00 

35.00 

5,00 

2000,00 

4000,00 

400,00 

0,00  0,00 

93.00  0,00 

18.00  0,00 

45.00 

45.00 

0,00 

5.00 
10,00 

1.00 

-4.00 

1,00 

1,00 

-45.89 
45,66 
• ••• 

0.84 

11.72 

2000,02 

4000,15 

«••• 

NUMBER 

or  states 

S 7776 

I 

N P U T ERROR 

S T A T 1 

: s t i c 

S 

SCOUT  NAVIGATION 
(METERS) 

X Y 

alt  heaojng 

(«)  (OEC) 

PITCH 

(OEG) 

roll 

(DEC) 

REL.  AZ. 
(RAO) 

REL,  EL. 
(RAD) 

RANGE 

(METERS) 

MEAN  0,0  0,0  0.0  0,0  0,0  0,0  0,0000  0,0000 

STANDARD  DEVIATION  55  ,g t 55,0.  ' 2*0  1,0  i.5 1,5  0,0002  0^000 2 


0.0 

10.0 


l 


computed  error  statistics  for  resulting  oensity  functions  about  target  x,y.z  coordinates 

x error  y error  z error 


~~  MEAN  ' 0,00  0.00  0,00 

STANDARD  DEVIATION  67,26  67,2.6..  26^87 
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